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Abstract 

We calculate the effect of zero-point fluctuations on the magnetostatic proper- 
ties of a quasi one-dimensional antiferromagnetic helimagnet with and without 
in-plane anisotropy. We find substantial reduction of the uniform suscepti- 
bility and the sublattice magnetization from their classical values and appre- 
ciable renormalization of the spin- flop field at T = 0. The magnetization 
curve varies nonlinearly at small fields and is described by a universal for- 
mula at high fields. The results are compared with numerical simulations on 
one-dimensional systems and measurements in the CsNiCla-type compounds. 
PACS numbers: 75.10Jm, 75.50Ee, 75.60.Ej 
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I. INTRODUCTION 



Recently, there has been growing interest in the behavior of one-dimensional and quasi 
one-dimensional antiferromagnets in strong magnetic fields. In the former case much at- 
tention was devoted to the destruction of the Haldane state of a spin-1 antiferromagnetic 
chain in applied field and to the magnetization process at higher fields [1-7]. In the latter 
case the interest lies mainly in properties of CsNiCls-type magnetic compounds which are 
good model Heisenberg systems. They contain a stacked triangular lattice of antiferromag- 
netically coupled 3d ions with the intrachain exchange constant J being several orders of 
magnitude larger than the interchain coupling J'. Because of such smallness of interchain 
interactions, quantum effects play an important role in the Neel phase. Extensive theoretical 
works have been presented to describe quantum renormalizations of the magnon spectrum 
at T = in these materials [9-11]. Though experiments indicate significant deviations in 
static properties of quasi one-dimensional antiferromagnets from predictions of the classical 
theory as well [12-18], only a few effects were considered theoretically [11,19,20]. 

Here we study the lowest order quantum corrections to the static characteristics of a quasi 
one-dimensional antiferromagnetic helimagnet: staggered magnetization, transverse and in- 
plane magnetic susceptibilities, spin-flop field, and full magnetization curve. We show that 
the magnetization process at high fields is the same for all low-dimensional systems, i.e., is 
independent from the nature of the ground state at H — 0. This fact is connected with the 
different corrections to the staggered and uniform magnetizations appearing in the spin-wave 
theory. 

The model is defined by the Heisenberg Hamiltonian for a system of N equivalent spins 
on a simple Bravais lattice with a single-ion anisotropy and the Zeeman term: 

n = Y: Jrj s, • s, -rDY^istf - 5:7H • s, , (1) 

i,i i i 

(in the following, we include the gyromagnetic ratio 7 = g/iB in the definition of H). With- 
out anisotropy and magnetic field the classical ground state of (1) is a plane spin helix 
Si — {S cos ^i, 5'sin ^j, 0), 9i — Q • r^, where the ordering wave vector Q corresponds to the 
minimum of the Fourier transform of exchange interactions: Jk = J2i J-i cxp(zk • rj) [21]. 

To obtain the first 1/5* correction to the classical energy of the ordered state we use the 
standard perturbation scheme in which spin operators are replaced by their bosonic equiva- 
lents by applying the reduced Holstein-Primakoff transformation — S — a'^a, — -\/2^a, 

2 



S~ — \/2Sa^ to each spin in its local coordinate system. The first order corrections ap- 
pear in this harmonic approximation due to the quadratic part of the boson Hamiltonian 
which corresponds to the quantized classical spin waves with zero-point oscillations. Diago- 
nalization of the boson Hamiltonian is trivial for a helimagnet without in-plane anisotropy 
(D > 0, H II z). Analytical expressions for quantum effects under these conditions are 
given in Sec. II. Less studied in the literature case of a bunched spiral in the presence of 
an in-plane anisotropy is considered in Sec. III. Treating perturbatively helix distortions 
we calculate the classical spin-wave spectrum and find the quantum effects for a spin helix 
with the easy- axis anisotropy and in-plane magnetic field. Comparison of the theory with 
experimental measurements in CsNiCla and CsMnBrs is presented in Sec. IV. 



II. AXIALLY SYMMETRIC CASE 



We first present standard spin-wave calculations for a helimagnet with the easy-plane 
type of anisotropy {D > 0) and magnetic field parallel to the symmetry axis z (see, 
e.g., Ref. [21]) and find analytical expressions for quantum effects in the quasi one- 
dimensional limit. For this axially symmetric case spins form an umbrella structure 
Sj = (5" cos q; cos ^i, 5" cos q; sin ^j, 5" sin a), where 9i — Q, ■ r^. By minimizing the classical en- 
ergy of the deformed spin helix one obtains an expression for the canting angle sin a — H/ Hs, 
which is vahd up to the saturation field Hg — 2(Jo — Jq + D)S. Above this critical field 
sin a = 1 and all spins are aligned parallel to H. 

The quadratic part of the Hamiltonian after the normal ordering and Fourier transfor- 
mation of boson operators is written in the form: 

n2 = ^NDS cos^ a + Y, [(^k + C'k)a,^ak + ^Bi,{ai^a_^ + a^a^)] , (2) 



with 



= -25 Jq + Scos^a{Jk + D) + + sin^a) (Jk+q + Jk-q) 



Bir — —S COS a 



Jk + D- |( Jk+Q + Jk-q) , C'k = 5" sin «( Jk+q - Jk-q) 



Note, that wc work in the extended zone scheme assigning wave vectors to the paramagnetic 
Brillouin zone. Coefficients Ak and i?k arc even functions of k, while Ck is an odd function. 
Performing the standard Bogolyubov transformation, which leaves odd-k terms unchanged, 
one can diagonalize the bilinear form of boson operators: 
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E h &k + 1) - 1^1. 



where the classical spin-wave frequencies are given by 

u;^ = JaI-BI + Q 



k • 



(3) 



(4) 



The easy- plane anisotropy creates gaps at k = ±Q, whereas spin waves at k — > remain to 
be the Goldstone modes up to the saturation field. 

We calculate the magnetization at T = taking the derivative of the ground state energy 
with respect to the magnetic field 



M 



H 



2( Jo -Jq + D) 



1 



2SN 



E 



Jo — Jq + D V ^k — -Dk 

The factor before brackets is the magnetization value in the classical theory which is renor- 
mahzed due to zero-point oscillations. The expression (5) corresponds to the field region 
< H < Hs, whereas in the completely aligned phase above quantum corrections disap- 
pear and magnetization saturates: Mg — NS. The average value of the local spin below Hg 
is given by the standard formula [21] 

1 _ / 



(5) 



1 - 



2SN 



E 



- 1 



*k 



(6) 



In the general case correction terms in (5) and (6) should be evaluated by numerical 
integration over the three-dimensional Brillouin zone. However, analytical expressions rep- 
resenting the leading logarithmic terms can be found in the case of quasi one-dimensional 
magnets. Let us calculate them for hexagonal CsNiCla-type antiferromagnets. Using the 
rectangular coordinate system, it is convenient to measure wave-vector components along 
the z and x axes in units of the reciprocal intersite distances c and a respectively, while 
the component in the y direction in units of 2/y/3a. Then, the Fourier transform of the 
exchange energy takes the form 



Jk = 2 J cos kz + 6 J'7k , 7k = (cos + 2 cos ^k^ cos k. 



yj ' 



(7) 



with each wave- vector component varying between and 2%. The helix vector of the trian- 
gular structure found by minimization of Jk is Q = (47r/3, 0, vr) in this coordinate frame. 

Keeping only the leading contributions in the quasi one- dimensional limit J ^ J', D, 
we omit terms proportional to J' and D in numerators of Eqs. (5) and (6) and otbtain the 
following expressions {Hs ~ 8JS): 



m = M/N 



H 
8J 



_ 1 r d-^i 
~2S J 



(— cos hr) x/l — COS /,•; 



(27r)^ ^1 + cos A;^ + pk - H^/{32J^S'^) cos A;, 



(8) 



'S'av — -S" 



1 /• d^A; 



1 -i7V(64J2'52)cosA;, 



1 — cos kz + Q'k) (1 + cos kz + Pk — H'^/ (327^5'^) cos kz 



(9) 



where 



3-^'/. o X D 3J' ^ 



The integration over kz in (8) can be done exphcitly for arbitrary pk and H. However, this 
would be an overestimation as we have aheady neglected linear and higher order terms in 
J'/ J. Therefore at H <^ Hg we integrate (8) keeping only logarithmic and zero order in pk 
terms. This yields: 

16J 



H r 1 
"^=8j1'-2.5 



{l - ^ [in ^ - 2 - (ln(l + 27k + D/3J' + HyASJfS'))] | , 



(10) 



where angular brackets denote averaging over the two-dimensional Brillouin zone. The 
sublattice magnetization is given with the same accuracy 0[{J'/ J) ln(J/J')] by 

5av = ^ {l - ^ [in 1^ - TT - i(ln(l - 7k)) - i(ln(l + 27k + D/3J' + HyA8JJ'S'))^ | . 

(11) 

The main part of quantum effects for both m and S^cv is represented at = by the 
large logarithm ln(16J/3J'), which is in agreement with a qualitative suggestion of Ref. [15]. 
Numerical constants yield in the isotropic case only small corrections: 

(ln(l + 27k)) = -0.452 , (ln(l - 7k)) = -0.176 . 



At low fields the magnetization is a linear function of H with the susceptibility per site 

16J 



_ 1 r 1 



{^-2^h3J^-2-(^^(l + 27k + W))]}. 



(12) 



The additional correction terms in (10) and (11), however, grow with increasing field and 
give nonlinear dependences of magnetization and spin length. This nonlinear increase of m 
and S'av means suppression of quantum fluctuations in magnetic field, which must completely 
disappear in the coUinear phase at H > Hg. 



In the isotropic case Eqs. (11) and (12) give good coincidence up to 2-3% with the results 
by Ohyama and Shiba [11], who estimated integrals in (5) and (6) numerically. Note, that 
all the derived expressions are also vahd for the easy-axis helimagnet at fields above the 
spin-flop field (H || z), when spins lie in the basal plane and the axial symmetry is restored. 
In this case the constant D in (10) and (11) has the negative sign as in the Hamiltonin (1). 

At high fields the magnitude of quantum corrections to m becomes independent of the 
interchain coupling. Crossover between two regimes takes place at a characteristic field 
He — y/A^J'JS"^. At He <ti H < Hg one can neglect the term pk in (8) and obtain after 
integration over the magnetization 



m 
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1 - 



where f(H) = < 




f{H) 

l\, H< HJV2 , 

HJV2 <H<H, 



(13) 



As for the staggered magnetization S'av, it remains finite in the field region < H < Hs only 
if J' ^ 0. 

According to (13) the magnetization approaches its saturation value with the asymp- 
totic behavior M — Ms[l — ^^2(1 — H/ Hs)/t:S]. As we neglected the interchain coupling, 
this result should be also vahd for a single antiferromagnetic chain. Hodgson and Parkin- 
son [3] using the Bethe ansatz approach found a different critical behavior for the chain 
magnetization: 



2^1 - H/HjnS 



(14) 



The reason of this discrepancy can be explained as follows. In the bosonic representation 
the Hamiltonian (1) is equivalent to a gas of Bose particles interacting via a (5-potential. 
Their interaction proportional to 1/S is neglected in the harmonic approximation. As mag- 
netic field approaches the saturation field, the boson density tends to zero. In the low-density 
limit the average kinetic energy of bosons is small and they behave as if the potential has 
a hard core. Therefore, in one dimension at fields sufficiently close to H^ the interaction 
between spin waves becomes non-negligible for any value of S and gives the modified asymp- 
totic behavior [6]. The width of the critical region of quantum fiuctuations decreases with 
increasing S. 
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We compare our expression for the magnetization (13) with the results of exact numerical 
diagonalization of the Hamiltonian for finite systems. As can be seen from Fig. 1, the formula 
(13) is in remarkable agreement with the magnetization curve for a spin-1 antiferromagnetic 
chain found numerically by Sakai and Takahashi [5] except for a low-field region where the 
nonmagnetic Haldane phase is developed and a close vicinity of Hs- Thus, we conclude that 
at the intermediate field region the magnetic field suppresses zero-point fiuctuations and 
increases the accuracy of the spin-wave theory. The magnetization curve for a system of 
coupled antiferromagnetic chains does not depend at high fields on the nature of the ground 
state (which is determined by J' and S) and is given by (13). However, in the narrow field 
region near Hg the interaction between spin-wave modes becomes important and changes 
the asymptotic behavior of M{H). The size of this region is small already for 5" = 1 and 
the harmonic approximation of the spin-wave theory gives good account for most of the 
experimental results at high fields [1,2,18]. 



III. IN-PLANE ANISOTROPY 



Consider further quantum effects for a Heisenberg helimagnet with an in-plane 
anisotropy. The easy-axis anisotropy or magnetic field applied in the spin plane perturb 
spin helix adding higher-order harmonics nQ [21,22]. Derivations of quantum corrections 
are similar for all cases. We present in detail calculations for < and H — listing in 
other cases only final expressions. 

If the anisotropy constant D in (1) is negative, spins in the spiral structure rotate in the 
vertical plane: Sj = (S'sin6'i, 0, S'cos^j). As the exchange interactions are much stronger 
than the relativistic anisotropy we can consider only the first anharmonic term 

^, = Q-r,-(^sin(2Q-rJ . (15) 

Making substitution of (15) into (1), we find the following expression for the deviation angle: 

^=\D\/{J,q-Jq) . 

The formula (15) is valid for both incommensurate and commensurate antiferromagnets (if 
only 3Q 7^ ±Q). In the latter case there should be an additional constant term in (15) which 
describes lock-in of a spin helix with respect to the crystal axes. However, we can neglect 
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it, because the commensuration energy appears only at the N/2 order of the expansion in 
\D\, where A'"Q = G, G being the reciprocal lattice vector [22]. 

Taking into account (15) the quadratic part of the Hamiltonian is written as 



k 

+ C'k(ak+Qak-Q + c.c.) + i^k(ak+Qa-k+Q + ak-Q«-k-Q + c-c) 



(16) 



where 



^k 
5k 
Ck 



-2JciS + S 



k-Qy 



-s 



D\S + jlfS Jk — ^{Jk+2Cl + Jk-2Q) 



D^^C,.- \D\S . 



As was shown by Elliot and Lange using general symmetry arguments [23] , a helimagnet 
with in-plane magnetic field always have a zero-frequency mode for spin rotations inside helix 
plane. Their theorem fulfills obviously for another types of in-plane anisotropy. We illustrate 
this by finding in the first order in D analytical expressions for the classical spectrum and 
the ground state energy of an easy-axis helimagnet. 

The Bogolyubov transformation can be used to diagonalize the first two terms in (16). 
After this procedure the boson part of the Hamiltonian (16) is reduced to 



k 

-k+Q + ^k-Q^-k-Q + C.C.) I , 



(17) 



where 



"^k 



Ai-Bl, tanh 2^k = -^k/Ak , 



Ck = Ck cosh('0k+Q + V'k-q) + -Dk sinh('0k+Q + V'k-q) , 
Dk = CkSinh(V'k+Q + V'k-Q) + -DkCOsh(V'k+Q + V'k-Q) ■ 
Note, that constants cu^, Ai^ and are given by the same expressions as in the case of 
the easy-plane anisotropy (2) but with Des = \D\/2. Attempt to find eigenfrequencies of 
(17) leads to an infinite system of coupled equations. However, the diagonalization may be 
simplified since the amplitudes Ck, -Dk are of the order of \D\S everywhere in the Brillouin 
zone except for vicinity of the point k = 0, where the Bogolyubov transformation becomes 
nonanalytical in the limit |D| 0. Therefore umklapp terms in the Hamiltonian (17) 



are important for the calculation of the spin- wave spectrum in the first order in \D\S only 
near k = ±Q, where they describe interaction of four nearly degenerate magnons with 
propagation vectors ±Q ± q. 

In the following calculations we use the method proposed by White [24]. We define 
the column vector of bosc operators Xq = 

constant term — |Z)k^k, the Hamiltonian (17) can be rewritten as 



q+Q' VQ'^-q+Q'^-q-q)- ^hcu, excluding the 



^-E^q^q^c 



(18) 



where 



q+Q 



Co 



V 



^q 





^q-Q 











q+Q 



(19) 



q-Q / 



Following the above arguments, we replace the eigenvalue problem for the infinite matrix H 
by diagonalization of 4 x 4 matrix (19). The transformation = Sy]^ that diagonalizes Hq 
must preserve the commutator = g, where the metric matrix g has only diagonal 

elements: (+1,-1-1,— 1,-1). Solving equation 



S ^gHqS = gQn , 

we obtain the diagonal matrix Qh with two pairs of equal eigenfrequencies 
(-^^r - ^[«q) + K-q)^] +Cl-Dl± {i[«Q)^ - (a;J_Q) 



(20) 



212 



q+Q 



(21) 



which can be mapped onto two different magnon branches in a half of the Brillouin zone. 
At k = Q the frequency gap of the first branch is loq — 25'[|D|(Jo -|- J2Q — 2Jq)]^/^, while 
cuq — 0. Thus, the classical spectrum of a helimagnet with the easy-axis anisotropy contains 
two quasi Goldstone modes at k = and Q, where the additional zero- gap mode results 
from the absence of the basal plane anisotropy in the Hamiltonian (1). Wc emphasize again 
that this conclusion is valid for a commensurate helimagnet only in the first order in \D\. 
For example, the second mode at k = Q in triangular antiferromagnets acquires a frequency 
gap proportional to [22]. 



The energy of zero-point fluctuations is determined after the last diagonalization by 
E^^^ — I Z^k;>o('^k +<^k)- However, the umklapp terms in (17) yield correction to the energy 
of the ground state only in the second order in anisotropy constant and, therefore, should 
be omitted. The ground state energy in the first order in |D| can be obtained after partial 
diagonalization: 

^g.s. = JqS{S + 1) - l\D\S{S + 1) + i E^k ■ (22) 

k 

To calculate the staggered magnetization with the same accuracy one has, in contrast, to 
determine explicitly the transformation matrix S. Expressing operators Oj via new boson 
modes, it is possible to find the average spin length including its modulation [~ cos(2Q • r^)] 
in a bunched spiral. These calculations cannot be done analytically and require numerical 
study which was presented in the case of triangular antiferromagnets by Watabe et al. [20]. 
Note, that for a small anisotropy (the case with which we are only concerned here) the 
sublattice magnetization can be roughly estimated as in the isotropic case by using Eq. (11) 
with D^Q. 

Expanding (22) up to the first order in we get the anisotropy energy of a triangular 
antiferromagnet renormalized by zero-point oscillations 



... \D\S{S+l) \D\S \^ 16 J ^ 
Ei"^ = + ^—^ In — - - (ln(l + 27k)) 



47r 



3J' 



(23) 



If magnetic field is applied along the easy axis, the vertical position of the helix plane becomes 
unstable at high fields. Magnetic energy exceeds the relativistic anisotropy, and above the 
spin-flop fleld H^i spins lie in the basal plane. In the classical picture, anisotropic part of the 
energy for this conflguration E'^ vanishes. However, zero-point motion of quantum spins 
makes E'^^ 0- Expanding the ground state energy (3) in powers of D, we have in the flrst 
order 



\D\S 
27r 



ln|^-(ln(l + 27k)) 



(24) 



The energy difference for vertical and horizontal configurations of polarization plane £^an — 
E^ — E^ will be used below for the calculation of the spin-flop fleld. 

Consider now the effect of magnetic fleld applied along the easy axis. The distortion of 
a spin helix is described in this case by 

= Q • r j - Q;sin(Q • rj - (^sin(2Q • rj 
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with a. — H / {Jq-\- J2Ci — 2Jq)S [21,22]. As we explained above, the spin harmonics generated 
by single-ion anisotropy, being important for calculations of the correct classical spin-wave 
spectrum, can be neglected for the estimate of the ground state energy in the first order in 
D. For the latter it is necessary to take into account only changes in constants and Sk. 
The new spin harmonic generated by magnetic field gives 

= (^Sk = -\a'^S [Jk-Q + Jk+Q - Jk - |(Jk-2Q + Jk+2Q)] ■ (25) 

Magnetic field also generates umklapp terms in the boson Hamiltonian with a momentum 
transfer ±Q. These terms are proportional to a and, therefore, should be taken into account 
for calculation of the ground state energy with the accuracy 0{a^) ~ H'^. Inspection of 
corresponding expressions reveals, however, that umklapp terms have an additional smallness 
of the order J' / J in the quasi one-dimensional limit. Having in mind this case, one can 
calculate magnetic contribution to the ground state energy without umklapp terms as in 
the above analysis of the easy-axis anisotropy. The uniform magnetization found in this 
approximation is 
H 



m 



2( Jo -|- J2Q — 2 Jq) 



^ _ 1 Jk + |(Jk+2Q + Jk-2Q) — Jk-Q — Jk+Q / Ak — -B; 



2SN Y >^o + J2Q - 2 Jq V ^k + 5k 

(26) 



1 r 1 

XII 



Analytical estimate of the integral in (26) yields the following expression for the in-plane 
susceptibility of a quasi one- dimensional triangular antiferromagnet 

16J 

<; I In 

16J 

Note, that quantum corrections break the relation x± = X\\y which holds for a classical 
antiferromagnetic helimagnet in the quasi one-dimensional limit. 

With the help of (12), (23), (24), and (27) we find the critical field of the spin-flop 
transition as a field at which anisotropy and magnetic energies of a quantum helimagnet are 
equal: 

^2 _ .^p .02 1 + 1/S- 3/(27r5)[ln(16J/3J0 - (ln(l + 27k))] 

° ■''^ 1 - l/(27r>5)[ln(16J/3J0 - 2 - 2(ln(l + 27k)) + (ln(l - 7k))] ' ^ ^ 

Let us list at the end the results for two another important cases, which can be derived 
in a similar way: (i) D < 0, H _L z, the magnetization is given by the formula for the 
easy-plane anisotropy (10) with the substitution D —>■ \D\/2; (ii) D > 0, H _L z, in-plane 
susceptibility in this case coincides with (27). 
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IV. DISCUSSION 



We compare our theoretical results with measurements on two isostructural quasi one- 
dimensional hexagonal antiferromagnets CsNiCls and CsMnBra, for which a lot of experi- 
mental data are available including the zero-point spin reductions [12,13] and precise magne- 
tization curves [14,15]. Both compounds are considered to be ideal Heisenberg systems de- 
scribed by the Hamiltonian (1) and differ only in sign and relative strength of the anisotropy 
constants. 

First we consider the case of CsMnBra, which has a large spin {S — 5/2) and the 
easy-plane type of anisotropy. The following values of the constants in the Hamiltonian 
(1) were obtained by fitting the measured magnon frequencies to the classical spin- wave 
dispersion [13,14]: J = 214 GHz, J' = 0.5 GHz, D = 1.95 GHz, and g ^ 2. It should be 
mentioned here that the magnon spectrum is also subjected to quantum renormalizations. 
If they are not too large the functional dependences of the magnon frequencies are expected 
to preserve the classical form, which involve, however, renormalized values J, J', and D. 
Major renormalizations in this case are found to be [25] 

For S — 5/2 these corrections are not so important and one can use the above set of constants 
to fit the experiment. With the help of (11) we obtain the average value of spin Sg^ — 
1.82, which is not far from the experimental value S'I^p = 1.65 [13]. However, the notable 
discrepancy about 10% is quite surprising. It cannot be explained by the renormalization of 
constants (29) and shows the limits of the harmonic approximation to describe a large spin 
reduction even for S = 5/2. Next we calculate the magnetization curve for H || z using the 
expression (10). The result is presented in Fig. 2 (solid line) together with the experimental 
data [14]. Theoretical curve coincides nicely with the experiment reproducing the unusually 
low value of the susceptibility and a small upturn curvature of the magnetization. Better 
agreement between the theory and experiment for the uniform magnetization than for the 
average spin value is clearly expected for the quasi one-dimensional antiferromagnet with a 
large easy-plane anisotropy D > J'. In the one-dimensional limit (J' = 0) the average spin 
value calculated in the linear spin-wave theory diverges showing its inapplicability, while 
the magnetization (10) remains to be finite. For H _L z we calculate only the slope of the 
magnetization curve using the expression (27) for the susceptibility. The resultant curve 
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(dashed line in Fig. 2) again coincides well with the experiment reproducing the quantum 
reduction of the susceptibility. At higher fields for this orientation of H the distortions of 
the triangular structure become important which requires to take into account the umklapp 
terms. At H > He the spin structure transforms into a coUinear one [14]. 

Next we consider CsNiCls which has the easy-axis type of anisotropy and much more 
pronounced quantum effects due to S* = 1. In order to calculate the magnetization for 
H ± z we use the expressions (10) with Dgs — \D\/2. For H || z the magnetization is 
given by (26) below the spin-fiop transition and by (10) but with negative D above iJgf- In 
the case of CsNiCls there is a serious problem to determine values of the constants in the 
Hamiltonian (1) from the spin-wave dispersions since the classical expression (22) does not 
reproduce the dependences observed on the experiment [8]. The discrepancy is small for 
the magnon dispersion along chains, and one can estimate from the zone boundary energy 
J = 345 GHz [8]. With this value of J the fit of the low-frequency magnetic resonance data 
and the spin-fiop field value to the classical expressions yield J' = 8 GHz and D = —0.6 GHz 
[26] . If we calculate the magnetization curves using this set of parameters, our results show 
rather poor agreement with the experimental data. Large part of the error comes from the 
use of the renormalized values for J, J', and D instead of the bare Hamiltonian constants. 
For S — 1 the renormalizations are quite important: according to (29) the bare value of J 
in this case is about 18% smaller than J, while D is twice larger than D. 

On the other hand, the bare value of the exchange constant J can be determined from 
the spin-fiip field Hg which is not renormalized due to the absence of quantum fiuctuations 
in the saturated phase at H > Hg. Recent high-field magnetization measurements [18] 
gave J — 275 GHz and g — 2.15. We used these values together with J' — % GHz and 
D — —1.6 GHz to calculate theoretical curves which are shown in Fig. 3. The choice of the 
bare anisotropy constant is rather arbitrary since it cannot be calculated from the spin-flop 
fleld: for such large value of the ratio J/ J' Eq. (28) fails due to the negative sign of the 
right hand side. Instead, we have chosen the value of D that gives reasonable description 
of the difference of magnetizations for the two perpendicular orientations of magnetic field 
at if > Hsi- Note, that this difference is connected with a zero-point contribution to the 
anisotropy energy, analogous to the case of CsMnBra above He [14,16]. Our calculations 
with the above constants presented in Fig. 3 give adequate values of the susceptibilities 
but distinct difference in the field dependence of the magnetization. We note that the 
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disagreement between the classical spin-flop field depicted in the Figure and the experimental 
value does not make any sense since both the classical VTODJS^ and the renormalized (28) 
expressions for Hgf are inapplicable for CsNiCls. This fact together with a poor agreement 
of the nonlinear growth of the magnetization point to the insufficiency of the harmonic 

approximation and need to develop a more accurate treatment including higher-order terms 
in 1/5*. Finally, the sublattice magnetization calculated neglecting the single-ion anisotropy 
is Say — 0.62 compared with the experimental value S^^^ — 0.52. 

V. CONCLUSIONS 

We calculated the leading contribution of zero-point fluctuations to the static properties 
of a quasi one-dimensional antiferromagnetic helimagnet. Use of the symmetry of the spin 
structure for I? > 0, H || z and perturbative treatment of the anisotropy for D < 
allowed us to develop a spin-wave theory with only one kind of boson operators in the 
paramagnetic Brillouin zone and to obtain simple analytic expressions for quantum effect. 
Note, that the traditional approach in the framework of the multisublattice model leads to 
large mathematical comphcations and does not yield analytical results [20] , though it should 
be applied in the case of a strong Ising-like anisotropy. 

Comparison with the experiment show that theory give good qualitative description of all 
specific features observed in the real quasi one-dimensional antiferromagnets. They include: 
strong reduction of the ordered moments and the susceptibility, nonlinear growth of the 
transverse magnetization and its anisotropy in the easy-axis case above Hg-f. Accurate 
quantitative comparison between the theory and experiment is difficult and give ambiguous 
results since the bare constants involved in the Hamiltonian (1) are not known. 

Zero-point fluctuations also affect in a different way different magnetic characteristics. 
The calculated correction to the spin-flop field makes its value imaginary in the first order 
in 1/5" for 5" = 1 what signals about importance of higher-order corrections. Instead, the 
results of the harmonic approximation give quite good description of the magnetization 
curve at high fields even for a single chain. Corresponding expression (13) fails only in the 
close neighborhood of the spin-flip fleld Hg where it should be corrected via account of the 
magnon interactions. 
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FIGURES 

FIG. 1. Magnetization curve for a spin-1 antiferromagnetic chain. Squares represent the 
numerical results of Ref. [5]. Solid line is the renormalized spin- wave theory prediction (13), 
dashed line is the magnetization asymptotics near the saturation field [3] . 

FIG. 2. Magnetization curves in CsMnBrs. Closed and open circles are experimental data 

per mole for H _L z and H || z, respectively [14]. Solid line is calculated using Eq. (10) with the 
parameters J = 214GHz, J' = 0.5GHz and D = 1.95GHz; dashed line is the low-field asymptotics 
with the same set of parameters. 

FIG. 3. Magnetization curves in CsNiCla. Closed and open circles are experimental data per 
mole for H _L z and H || z, respectively [15]. Solid lines are the renormalized spin-wave theory 
predictions with the parameters J = 275GHz, J' = 8GHz and D = — 1.6GHz. 
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